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CANONICAL GENERAL RELATIVITY: 

MATTER FIELDS IN A GENERAL LINEAR FRAME 

M. A. CLAYTON 

Abstract. Building on the results of previous work j^, we demonstrate how matter fields are 
incorporated into the general linear frame approach to general relativity. When considering the 
Maxwell one-form field, we find that the system that leads naturally to canonical vierbein gen¬ 
eral relativity has the extrinsic curvature of the Cauchy surface represented by gravitational as 
well as non-gravitational degrees of freedom. Nevertheless the metric compatibility conditions are 
undisturbed, and this apparent derivative-coupling is seen to be an effect of working with (possibly 
orthonormal) linear frames. The formalism is adapted to consider a Dirac Fermion, where we find 
that a milder form of this apparent derivative-coupling appears. 


I. Introduction 

This work is essentially an extension of that presented in |^], showing that the Hamiltonian 
description of matter fields minimally coupled to general relativity (GR) is compatible with the 
general linear frame formalism. In |5[ we gave a detailed description of the surface geometry and 
canonical structure of general relativity in a general linear frame, using two different choices of 
the diffeomorphism constraints: the ‘unprimed’ constraints which are compatible with a coordinate 
frame gauge choice, and the ‘primed’ constraints which are compatible with an orthonormal (or 
Lorentz) frame gauge choice. 

Briefly, we began with an action that allowed one to treat the metric and vierbein fields inde¬ 
pendently. Then, specializing to a surface-normal frame, we developed a Hamiltonian formalism 
that allowed one to consider the spatial frame independently from the spatial metric degrees of 
freedom as initial data, each with its own conjugate momenta. (Note that we have followed Q and 
specialized somewhat by choosing a surface-normal frame and assuming the existence of a spatial 
coordinate frame, choosing a minimal representation of the evolution of the spatial hypersurface 
in spacetime—for alternative points of view, see [^, |T^. Operationally this should be seen as no 
particular limitation since any numerical scheme would certainly begin from such a point of view.) 
We naturally end up with atlas fields (A, A“, that enforce the Hamiltonian, momentum and 
frame constraints respectively. The latter constraints generate infinitesimal frame transforma¬ 
tions, and are common to the two choices of diffeomorphism generators considered: the unprimed 
T-L and THul primed li! and 2i.'a constraints, related by nontrivial factors involving 

We shall see that, as expected, the addition of matter fields poses no great difficulty. The some¬ 
what surprising feature that appears is the breakdown of DeWitt’s “Riemannian Structure” 0 , as 
matter fields enter into the definition of the momentum conjugate to the frame fields (this should 
appear for any tensor field that behaves non-trivially under frame transformations). This is opera¬ 
tionally due to the presence of derivatives of the frame field in the matter action, entering through 
the non-vanishing structure ‘constants’ CabP which appear. Note that even though the gravita¬ 
tional and matter sectors are mixed, we show that the evolution of the spatial metric is identical 
to that required by the compatibility equations (i.e., the definition of the extrinsic curvature is 
unaltered) and the field equations are equivalent to those derived in a coordinate frame. Thus this 
form of derivative-coupling is benign in the sense that the light cone structure of the theory is 
undisturbed 0 — trivially because it is a coordinate frame theory that has been re-written. De¬ 
spite naive expectation, this works because the matter stress-energy tensor remains independent of 
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second-derivatives of the gravitational variables. We also find that the Einstein-Dirac (E-D) system 
is also derivative-coupled in this sense, and may be easily treated via a restriction of the formalism. 

The “unprimed” system of constraints is closely related to the coordinate frame approaches: 
evolution and Lie transport act solely on the components of tensors and the frame plays a passive 
role. To be specihc, infinitesimal spatial diffeomorphisms are dehned to leave the frame unaffected 
and change the tensor components in the usual manner [^] ^ : Ta —> £^\T]a. Writing Einstein’s 
equations in terms of the components of the metric 7 ^ 6 , frame E^a and extrinsic curvature Kab of 
a spacelike hypersurface S, we have 


d±[7]afe = 2iLafe, d^[E*a]=0, (la) 

<ll.[K]ab = -Rab + VaVb[N]/N - KKab + 2KacK% + (lb) 

where d^ is the general linear frame generalization (see the discussion in Section III-B of Q for more 
details) of the surface-covariant normal derivative operator |l^. We have added a matter stress- 
energy contribution Tab '■= Tab — ^EabT to for later reference, and we fix IGttG = c = 1. 
Lowercase Greek and Latin letters represent spacetime and spatial coordinate components and 
upper and lower-case Roman letters indicate spacetime and spatial frame components respectively. 
We also use a for the spacetime covariant derivative and Va for the intrinsic surface covariant 
derivative; more details may be found in 

In the Hamiltonian system the action of £ is represented on phase space by the momentum 
constraint as T^[-] —> We use the notation, for example, •= 

and {•,•} are the Poisson brackets given by (32) of [^. The Hamiltonian constraint represents 
d_L on phase space, and reproduces (||) weakly: 

{iab,n^V]} = VKb, = 0 , (2a) 

[kab^n^^U]] = -fRab + VaVb[/] - fkkab + 2fkack^ " (2b) 

{Kb,2i^^[f]} = -fRab + V,Vfe[/] - fkkU + (2c) 


(We denote the symmetric and antisymmetric projection on any pair of indices by, for example, 
T{ab) •= \{Tab + Tfta) and Ty^b] '■= \{Tab — Tba), densities with respect to the spatial metric by 
boldfaced symbols, and densities with respect to the determinant of the spatial vierbein with an 
underline.) We found that the extrinsic curvature Kab was equivalently represented on phase space 
by the two tensors (vr^^ and are the momenta conjugate to 7 “^ and respectively) 

kab ■— ~i'^ab ~ jTafeT), kab b)i (3) 

and the weak equivalence of (^) and (|^ guarantee that the evolution equations reproduce Ein¬ 
stein’s equations (|l]). Since we are working with a general linear frame on S we also have the gener¬ 
ators of frame transformations which act, for example, as |Tn, ^ [w'^fe]} = Ai,[T]a = -LV^aTb, 
and satisfy the Lie algebra of gl(3,M) 


Note that = —2{kab — 'Jabk) + 2{k'^^^ — 'jabk'), which guarantees that kab ~ The constraint 
algebra consists of (Q) and 

= 0, {2ia[n,J\[^’’a]} = X T A^[n]a, (5a) 

{mfimg]} = (/Vab] - 5V„[/])7“'2i6, (5b) 

{mflUaid"]} = [ dxf_£^[Ul {'HanUbW]] = [ d^xr£^[n]a, (5c) 

Js Js 

which was also derived in |^] using the geometric arguments of Teitelboim [l^ , and which is a 
fairly straightforward generalization of that found elsewhere (see for example [13(| ). 

We then noted from the form of (|T^) or (|^ ) that these generators are not very convenient if one 
wants to consider the restriction to orthonormal frames on S, since if one has chosen 'jab = dab, the 
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actions of d_L and £ (and therefore !£ do not leave this form invariant. There are, however, 

alternative representations of infinitesimal diffeomorphisms and normal evolution in which the 
spatial frame plays a more active role. In particular, we represent an infinitesimal diffeomorphism 
by £' which acts on a frame as £'^E^a\ = and on the components of tensors as the 

covariant derivative Th[T]a = so that the action on the tensor itself is identical to that of 

£: £^[Ta9°'] = Th[Ta0“]. The result is that the action of £' on the components of the spatial metric 
vanishes due to metric compatibility, and so its action is consistent with the limit to orthonormal 
spatial frames—note though that it is no longer consistent with the limit to a coordinate frame. 
Similarly we define the operator d(j_ to act of the frame as d(j_[£'*a] = Aj^[E^a] = —K^aE^b and on 

the components of tensors as d(|_ = d_L + A^. Using these, Einstein’s equations (||) appear as 

dl [74 = 0, d'^[E\]=Aj^E% = -K\E%, (6a) 

d±[KU = -Rab + VaVb[A]/A - KKab + df^b- (6b) 


The action of d^^ and £' are represented on phase space by the primed Hamiltonian constraint 
EC and the primed momentum constraint EL '^ respectively (which correspond to a particular choice 
of mixing of the unprimed constraints with the frame rotation generators), giving 

[lab.'lH^^lf]] = 0, = -k'\E\, 

= -fRab + VaVfe[/] - fk'kab - 

{k'ab,2i'^V]} = -fRab + VaVbif] - fk'k',, - kf^abU'^^ - ' 

The constraint algebra for this primed system consists of (® combined with 


(7a) 

(7b) 

(7c) 


{2i![f],J\[co’^a]} = 0 , {2iLW%J\[co\]} = l^d^xN'^Acbma, 

{n’[f],n'[g]} = [ d"x(/v,b]-5V„[/])(7“'2i'6-i?Vb[j'][“']), 

JT, 

{THainiUbia’’]} = - jj^xrg^R%abJ\, 


(8a) 

(8b) 

(8c) 

(8d) 


which was also derived in |p. 

As we shall see, including matter presents no great difficulties. In Section |n| we consider a scalar 
field which, since it is chosen to transform trivially under frame rotations, is a brief demonstration 
of how matter fields fit into the generalized structure. In Section HI, a one-form (gauge) field is 
introduced, which transforms non-trivially under spatial frame rotations. In addition, a constraint 
which generates U(l) transformations on the Cauchy surface appears, resulting in an extended 
diffeomorphism algebra. Here we find a nontrivial mixing of the Maxwell canonical pair with that 
of the spatial vierbein—mixing that, at least superficially, resembles derivative-coupling. We will 
see that this mixing is merely a manifestation of the gauge choice and is not “true” derivative- 
coupling since the metric compatibility conditions are not altered. Finally, in Section IV, the 
spatial frame is constrained to be orthogonal and a Dirac field is introduced. 

In all three cases the matter action is added to that of vacuum GR: and since 

both sets of vacuum GR constraints {EC^^^, Ei^^) and ki'a^) are equivalent to {—2G^±, —2G^a) 

{Gab is the Einstein tensor), we shall see that the additively combined vacuum GR and matter 
constraints are equivalent to adding (R-*-^, R-’-q) to these. The matter stress-energy tensor is con¬ 
veniently derived from the matter action via 


Tab 


2 f 6S^ 1 „ „CD \ 

E “ 2^ABg j 


1 

7 ^ 


^AcE^B 


6S^ 


gABE^C 


6S^ \ 


(9) 
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the equivalence of these two forms is a basic consistency requirement for the generalized treatment 
considered herein. By an argument of Floreanini and Percacci we know that once the matter 
action has been properly written in a general linear frame, these forms will be consistent. Note 
however that in (P) we are treating the metric and vierbein degrees of freedom in the matter action 
independently, which will not work with the E-D system where only the second form is applicable. 
We will show (briefly) that the matter field equations and the stress-energy tensor determined from 
both forms of (|^ are properly generated via a variational principle, and in the canonical picture 
that the action of the constraints on phase space properly generates spatial diffeomorphisms once 
matter has been incorporated into the canonical formalism, and (§) and (|^ are properly extended 
to reproduce (||) and i) respectively with the appropriate stress-energy tensor. 

One of the unfortunate aspects and at the same time strengths of diffeomorphism invariant 
theories is the great arbitrariness in parameterizing configuration and momentum space fields; any 
“all-encompassing” formalism would have to include canonical transformations that relate different 
parameterizations. The content of the ‘geometric’ derivations of the constraint algebras 11,11, i 
is that the algebra is fixed by the reduction of the spacetime equations to the evolution of spatial 
quantities; that is, the algebra is a reflection of the spacetime diffeomorphism invariance. The 
details of the algebra will change depending on the parameterization of phase space and the choice 
of constraints, however are always of a fixed form once these choices are made. Thus the point of 
this work is not to advocate any particular parameterization, merely to make it clear that there 
is actually a wider range of possibilities than heretofore considered. Within this wider class one 
finds that the Einstein-Dirac system may be treated in a particular limit, and not as an awkward 
construction that seemingly must be invented solely to deal with Eermion fields. Thus we find 
that the Einstein-Dirac (E-D) system is operationally no different than the Einstein-Klein Gordon 
(E-KG) or Einstein-Maxwell (E-M) systems; there is no more derivative-coupling in the E-D system 
than in the E-M system. Indeed the mixing of canonical variables is necessary in order to retrieve the 
correct surface frame transformation generators, while retaining the metric compatibility conditions 
at the level of the field equations. 


II. Scalar Eields 


The scalar field example that we consider in this section is a trivial extension of the vacuum GR 
results, included for completeness and to give a straightforward example of how the extension of 
the results of Q outlined in Section | to “GR-|-matter” systems proceeds. Its simplicity is due to 
the fact that we have chosen canonical coordinates that do not transform under changes of spatial 
frame. Note that this is not necessary; we could just as well chosen to parameterize the scalar 
held action in terms of a densitized scalar held, in which case some of the structure that we will 


encounter in Section ID would have appeared. 

The scalar held Lagrangian density with self-coupling potential V[(p] is 

£0 = - EM) 


from which we derive (from either form of 


^-^NY’’Vambm-NV[<P], 


1 


TIj, = Vb[</)] - -gABi 


CDA^ 


Vc[</)]VD[</)]+gABEM, 


T^^ = ^VA[4>]^VB[4>]-gABV[4>]. 


Using the held equations for (p 




= Eg^^^VA[^VB[4>]] -Ed^V], 


( 10 ) 

(lla) 

(llb) 

(12) 


it is straightforward to demonstrate that the conservation law b[T '^]^= 0 is satished. 
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From (1C) we see that the phase space of the scalar field is parameterized by a conjugate pair 
{(p,P) where the conjugate momentum is 


E 


P:= - N<^ea[cl)]) = Ee^lct)]. 


(13) 


As a result, the Poisson brackets are written as {F, G} = {F, G}gr + {-F, where {., .}gr are the 
GR sector Poisson brackets given in p] and 


{F,G}^ ■= 


6F 6G 


5G 5F 


S(j){x) (iP(x) 54>{x) dP_{x) 


(14) 


GRd 


is the additional contribution from the scalar field. 

Since neither the extrinsic curvature nor time derivatives of the spatial metric or frame appear 
in the Lagrangian density (IR), we see that the definitions of kab and in (^) are undisturbed 

Thus we 

(15a) 
(15b) 


b- 


and the frame rotation generators are identical to those in vacuum GR: ^ 

find that 

= \E{P^f + \E-i^^\Ia[(t>Wb[<P] + EV[cl)] = 

= = Z0V,[(/)]=T^\, 

are added to the constraints of GR. 

Hamilton’s equations for the scalar field 

{ct),H} = NP^ + N^Va[(^], {P4„H} = EVa[NY’’^h[(t>]] - N6^[V] + EVaW^P^], ( 16 ) 

are equivalent to (|^), and noting that {'7°'^,2^^[f]} = [/“]} = 0 and = 

[f°]| = 0, we see that neither of (^) nor ( [7a| ) are altered. Using (with similar definitions 
with primes for Poisson brackets of the primed constraints) 


= Y\[f]E\ - lY[f]E%, 

{lLab.U^[f]]=XM-labX[fl 

(17a) 

{p\,nt[f]} = w\[f]E% - 

{Aab.'ii^[n]=^b[f\-iabm, 

(17b) 

where 



w\[f] = w'\[f] = -rnf, 

Zab[f\ = Z'ab[f\ = 0, 

(18a) 

Xabif] = X^^,[f] = = - 

-\f{fab-labf)-\flahnt 

(18b) 

Yabif] = Y'M = -fTi = - 

f{fab - labf) - flabPt 

(18c) 

it is straightforward to show that the evolution of the extrinsic curvature is correctly gener¬ 
ated. That is, we find that and is appended to the right hand sides 

of (2b), (1^), (^), and (^ in accordance with (p]) and (|^. The unprimed and primed constraint 


algebras are identical to (44) and (49) respectively of with TC = and Tia = + 'kia', 

or, equivalently, (^0|) and (|4^) combined with = 0. 

III. Maxwell Fields 

The Einstein-Maxwell system is more involved. As the vector potential is a one-form field, 
it is associated to GLM through a (co-)vector representation of GL(4,M), and therefore results 
in a nontrivial contribution to J^'b- This means that the primed and unprimed cases may no 
longer be considered together; the unprimed system is considered in Section tll-A and the primed 
in Section III-B| . We shall see that although computationally involved, the results are a fairly 
straightforward extension of those of vacuum GR. 

The components of the field strength tensor F := dA = ^Fab^^ A 9^ are derived from the 
one-form A = AaO^ by 

Fab = b[^a = gaIAb] — eR[Ayi] — Cab^Ac, 


(19) 
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and the standard Lagrangian density is 

£A ^ -lEg^^g^^FABFcD = - iFabF^^) , (20) 

where Fab = '^a^b - ^b^a and 

F^a = {dt[Aa] - AbE\dt[E\] - £j^[A]a)/N - - A±Va[ln(7V)]. (21) 

The Maxwell field equations derived from (|^ ) are 

rcM 

-— = E^VB[Ff^ = Q. ( 22 ) 

OAa 

Making use of the variation of the field strength tensor with respect to the spatial frame: 

6eFab = 2^Vc[A\BE^^6E^^A\ + 2Ac^V yslE^ ^5E>^ a]1 (23) 

we find that the stress-energy tensor as determined by either form of (|^ is 

TXb = Tab = -g^^EAcEBD + \gABE^^FBc- (24) 


From the field equations (|^) and using the Bianchi identity a[F]bc+'^'^c[F]ab + '^'^b[F]cA = 0, 
we find that the conservation laws = 0 are satisfied. Note that despite the presence 

of the structure constants Cab^ in the field strength tensor (W), the variation of the action (^) 
using (^) does not lead to any terms containing the second-derivative of the frame. 

are determined as 


We see from (^) and {21) that the momenta conjugate to Aa 

6L^ 


= Y^F^b =■■ EX 


5dt[Aa] ^ 


(25) 


Therefore the phase space of the Maxwell field is parameterized by the conjugate pair {Aa,E)i nnd 
the Poisson brackets are extended by (summation over a is assumed) 



6F 6G 

SAa{x) 5P\x) 


6G 6F \ 

SAa{x) 5P“(x) ) • 


(26) 


However we also find that the definition of the momenta conjugate to EX a as given in equation (30) 
of II is extended to 

:= -2K\E\ - P^AbE\ + \A^PbE\, (27) 

which is the “geometrodynamical momentum” conjugate to the frame, the first term of which is 
the “gravitational momentum” |^. Thus we find that this definition mixes the gravitational and 
non-gravitational fields, and DeWitt’s super-metric is no longer block-diagonal. By defining 

Kafe := -\{P{aAb) - llabE^Aa), (28) 

we now have (c./. the discussion following (33) of ||) 

Kab ~ dkab T (1 F){kab T f^ab)j (29) 


where the unprimed case corresponds to a = 1 and the primed case to a = 0. In addition, the 
generators of frame transformations become 

JXb = ‘^YYbc - IL^t - P\E\ + p5t - P-Ab, (30) 

which correctly generate frame transformations on the vector field sector 

{Aa,j\[Yb]} = A^[7l]a, {E\J\[YX\} = A^[P]“, (31) 

and satisfy (H). 

The Lagrange multiplier A_\_ enforces the U(l) constraint: 

UXl) = -EValPW ( 32 ) 

and we can immediately compute 

{2i(j(i)[«],2i(j(i)[/3]} = 0, = 0, (33) 
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and the non-vanishing brackets 

=Va[«], {P“,2iu(i)[«]} =0, (34a) 

{p“D2iu(i)H} = - ^P^E\)Vb[a]. (34b) 

This additional constraint is appended to the standard forms of the Hamiltonian via an additional 
atlas field a := NA_i_ as (and similarly for the primed system) 


m = j d^x (iv?i+ iV“2ia + ^“6^^ + «2iu(i))- 


(35) 


To determine the form of the other constraints appearing in (^^ we follow the development in Q: 
choosing Kab = kah {a = 1) leads to the unprimed system that closely resembles the coordinate 
frame approach (Section [1II-A|) , whereas choosing Kab = k'ab + ~ leads to the primed case 

(Section III-B). 


Ill-A. The Unprimed System. Since for the unprimed system we have Kab = kab the results 
are straightforward to derive, merely pulling off the Hamiltonian and momentum constraints from 
the Maxwell Hamiltonian, resulting in the following additions to the unprimed constraints 


2^ : = kEP^Pa + lEF^’^Fab = 
m ■ = P^Fab - A^blPf = T^^a + ^aKu(l). 


From these it is straightforward to determine 

{Aa,nV]] = fPa, {P^,mf]} = -EVbifF^ 

{Aa,m[f]} = £f[A]a, {p^mif]} = E£f[pr, 

and using (|34|), Hamilton’s equations for the Maxwell field are found to be 

{Aa,H} =NPa + £j^[A]a + Af,[A]a + Va[a], 

{P“, H} = -EVb[NF^^] + E£j^[P]^ + A^[P]“, 


which, when combined combined with the U(l) constraint (32), are equivalent to (22). 
Using the definitions as given in (p^, we find 


w\[f] = -T/AfeP^] - rul 

Y\[f] = -ff^b - fPaPb - AbVaifF^^], 


Fab = 0, 

Xab[f] = -Iff^b- 


(36a) 

(36b) 


(37a) 

(37b) 


(38a) 

(38b) 


(39a) 

(39b) 


Noting that T := 'y°‘^Tab = it is straightforward to show that {kab,2i} = ^fPab — ^'Jab'F^, and 
therefore the correct extension of ( p1^ is generated. In order to check the extension of (^^, note 
that we now must compute + Kab,2i.[f]} which correctly results in the additional contributions 

yfab - IflabH^ - fkKab + 2fk(a‘'>^cb) + fk(a‘'J{^ab)V 

It is now possible to compute the algebra of the unprimed constraints, finding (Q), (|^), (|^, (^) 
and 


J s 

[mfimg]] = UXa[g] - gXa[f])i^\Ub - A2iu(i)). 


(40a) 

(40b) 


Note that the presence of the U(l) constraint is a new feature of the algebra not anticipated by the 
analysis of Q. 
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III-B. The Primed System. The primed constraints are 

^ ^ \EF^^Fab + {k\k\ - K^) + 2{K\k'\ - k'n), 

2i'a = P^^a[A\b, 


(41a) 

(41b) 


where note that they are related to (^) via combinations of the frame rotation generators as 
given in (45) of [||, and therefore the combined constraints are equivalent to the Einstein equations 


G'^A = A- Their action on the Maxwell phase space is found to be 


2 


[Aa^lH^f]] = fPa + = -EV^ifF^^] + /A^,^-[P]“, (42a) 

{^a,2i'a[r]} = f^b[A]a = £'f[AU {P^lUtin} = EVkif’^P^] = PflRV- (42b) 

Note that the actions (|^ and (|4^ are those of d_L and £ and d(j_ and £' respectively. Hamilton’s 
equations for the Maxwell sector are straightforward to find, and, as in the case of vacuum GR 
(see (52) of i), are weakly equivalent to (^) when one makes the replacements N = N', Na = N^, 
N\ = N'\ + Vb[NT + N'k\. 

Because we find {E^ai2i!^[f]} = the action of the Hamiltonian constraint on the frame 

becomes 


{E\,2^[f]} = fAy^-[E^J, 


(43) 


properly reproducing the evolution of the vierbein in accordance with (^). Again making use of 
the definitions (^) we find 


WU[f] = -faWt + V,[rP^bAa]] + V,[P'^A(6A)] - Ve[A7(6Pa)], 

Kbif] = -hfTab + Iflabib^'^bK^ - + flab{k'\K\ - k' k) 

yLbif] = -fTab - fPaPb + AbVcifFa^] 

- 2f{Ka''k',,^ - Kk'^^) - /7afe(K“fcK\ - K^) + f{Ka^ - 

Z'ab\f] = \{Vc[P^fiaAb)\ + V,[rPyAy] - V ,[A<^ f yPy]) , 

and noting that i^(a^J\^b)\ ~ \P^Pc.AaAb - \A^AcPaPb-, it is straightforward to show that the 
appropriate additions to (|fb|) and (Hi are generated from these. 

The primed constraint algebra consists of (|i, (|^ , (|8a|) , (|8d|) and 


(44a) 

(44b) 

(44c) 

(44d) 


{2iu(i)[a] 

,mf]] = 0, [mjiyWi'Kain] = - jj^xa£'/^U'hy)l 

(45a) 

[mfimg]} = 

^ {fVa[g] - gVa[f]){y'^^n\ - 

(45b) 

{mfiiHam = 

f d?X (/%[?i'] - [K']a + 2g^V,[f{k'ab + Kab)]^^^ 

J s 

^1). (45c) 


IV. Dirac Spinors 


In order to consider the introduction of a self-gravitating Dirac spinor, we need to do some 
preparatory work. The conditions for the existence of a spin structure on a manifold are known |9|, Q , 


however since we are primarily interested in the initial-value formalism, the result that any globally 
hyperbolic spacetime or spacetime that admits a Cauchy surface admits a spin structure |]T^ is the 
most relevant. The general theory of spinors in curved spacetime (see [^|) requires the introduction 
of an SL(2, C) (the universal covering group of the Lorentz group) principle bundle, and a Dirac 
spinor is associated to it through the product of a vector and conjugate vector representations of 
SL(2,C). Generally one introduces spin frames, and we would then need to develop formalism to 
deal with the initial-value problem [p^j. 
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Here we will pursue a more straightforward, operationally-oriented path. We note that a Dirac 
spinor is associated directly to the Lorentz group through a spinor representation |]^], and so we 
will reduce the general linear frame bundle GLM to the Lorentz frame bundle LM, with structure 
group 0(1, 3). (In fact, for the initial value problem we will be considering time and space oriented 
frames and so have in effect reduced this further to M with structure group the proper Lorentz 
group £i^(l,3) = {A G 0(1,3)1 det(A) = -|-l,A*^o > 0}.) Elements of GL(4,M) are specialized to 

0(1,3) as A^b and Ab^ = V^'^Vbd-^^c, where A^c^b^ — ~ B^c^- 

We assume that given a Lorentz transformation A, the spinor transforms as 


i; - S{A)i,, 


(46) 


the gamma matrices that satisfy the usual Dirac relations (we make use of the definitions and 
identities in Appendix 2A of |2^) 

{ 7 ^, 7 ''} = 2r?^^, (47) 

and satisfy A^bI^ = 5 '“^(A) 7 "^S'(A). Writing an infinitesimal Lorentz transformation as A^b ~ 
(5^ -|- uo^Bi ^B^ ~ “ ^^^B, and lo^b are the antisymmetric generators of so(l, 3), for the spinor 

representation we write S{A) « 1 — where oj := oj^bo'^a and := satisfy the 

commutators 


^CD] ^ _2i{^jAC^BD _ ^AD^BC ^ ^BD^AC _ ^BC^ADy 
We therefore define the action of the generator of frame rotations A on the spinor as 

(49) 

and using (48) it is straightforward to show that these generators satisfy the algebra 

(50) 

The covariant derivative operates on these spinors as 

^Va[V’] = &aW\ + \i^AC(^^Bi’ = caM + if^f’AV’, (51) 

and we note that [^Va, 7 ^] = 0 . 

The curved spacetime Dirac action is (see Section 7.10.2 of IITtII ) 

5° = J d'^x A — (52) 


where we employ the standard notation aV’ •= iA^Va[iA] — a['4’]'4’■ Variation of this with 

respect to and 'll: leads to the Dirac equation 


aW\ — "rmp = 0 , 


(53) 


and its adjoint respectively. As a consequence, the vector := ■ 07"^'0 constructed from a solution 
to the Dirac equation is covariantly conserved: a[J]"^ = 0, which in the surface-adapted frame 

becomes dtlEij:^if:] — EVa[N°"ij:^ip] + E'Va[Nip^a°‘ip] = 0 and has the usual quantum-mechanical 
probabilistic interpretation. Unlike the scalar and Maxwell field examples considered previously, 
the spinor is not dimensionless after one has chosen IGvrG = c = 1. Indeed, this length scale may 
be interpreted as the Planck length Lp, or h . Here we may either consider ip to have a dimension 
of length~^^‘^ or to have taken h = 1. Note also that in this section we will introduce neither 
7 ^ := E^aA^ nor the coordinate components of the spacetime metric. 

Using the variation of the ‘spin-connection’ with respect to the vierbein 

5f A = SE^aE^^^Fd 

+ [vbd'^VaIE^^SE^^c] - vbd^^c[E^^6E^a\ + vad^^b[E^^SE^c])(t^^, ^ 

and the second form of d), we find the stress-energy tensor 
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Once again we find that there are no terms containing second derivatives of either the spatial metric 
or frame, albeit trivially because the Dirac action only contains terms that are at most linear 
in first-order derivatives of these fields. The conservation laws a ~ follow from the use 

of (|^) , the Bianchi identities of the Riemann tensor and the wave equation satisfied by the Fermion 
field: B + m^ — \R)'il} = 0 , where we have used that ^Vs]V’ = \iR^ dab^^^ c'4^ = 

jiRABi^ and Rab<^^^ = —2R. 


IV-A. The Initial-Value Formalism. To specialize the primed system of vacuum GR to consider 
orthonormal frames is a simple matter of replacing jab by 6ab in all results and ignoring anything 
related to the ilLab) Part of phase space (ending up with a description of canonical tetrad GR 
similar to that discussed in Q). In addition, it will be useful to introduce (using the conventions 
of Section 1.1 of [22|) (5 = and = 7 *^ 7 “, in terms of which we have cr^“ = and the spin 
matrices cr“^ := -i-aloRiof’] satisfying the reduced version of 


cr , a 


We define the reduction of the operator 


+ + + (56) 

that generates infinitesimal 0(3) frame rotations by 


= ijOJil), where Cj = Wafecr 


ab 


(57) 


and the surface-covariant derivative as 


Va['!/^] = Cailp] + (58) 

where Ta := dbdX'ac^'^^- commutator of these derivatives acting on a spinor results in [Vq, Vfe ]'0 = 

liRab^ = \iRcdab(y‘'‘^'4’- 

The decomposition of the spin connection is given by Ua = i\Ta + and = 

^aacR" — i\C±_ {C±_ := 5bcC^a^a°‘^ and Cacl^ is given in equation (18b) of and it’s role in the 
surface covariant normal derivative is discussed in Section III-B of the same reference), so that we 
have 


ValV’] = V.fV'] + (59a) 

"VxlV’] = + \aaa^'>P + ijJEa^^t[E\]a^^^l,. (59b) 

Using these and introducing the “half-densitized” spinor fields (since in this section we have 
det( 7 a 6 ) = det((5afe) = 1, densities of this type will not appear and there should be no confusion 
with boldfaced quantities that appear elsewhere) 

■0 := 0^ := F^/^0^, (60) 


from (52) it is straightforward to deduce the curved spacetime Dirac Lagrangian 


L° = 


J d^x (^^0^ d t'4’ + f^iVF0'^Q;“ V a0 — mN'il^'^Prp 

[iV]60t-“'’ 


iFV“0^ V a0 + iV„ 


V - dFai5t[F;*fe]0V“'’0). (61) 

In order to perform the Legendre transform to obtain the E-D Hamiltonian, we need to choose 
an appropriate set of canonical coordinates in the Fermionic sector. In a flat spacetime this is 
straightforward: one finds that (discarding a total time derivative from the action) the real and 
imaginary parts of 0 are canonically conjugate, and one can introduce complex coordinates on 
phase space and use the quantum-mechanical symplectic form |]T] for a four component vector. In 
a curved spacetime one hnds that it is the densitized imaginary part of 0 that is conjugate to the 
real part, however introducing the half-densitized spinors ( |^ allows the flat spacetime construction 
to proceed in an identical manner. Using the inner product on S (written first in covariant form, 
and then specialized to a surface-normal frame) 


(0,0)s= / dJ:{x)nAip'y^(p = 


/ dxip^cf) 

Jt, 


(62) 
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to define the pairing and therefore the quantum-mechanical symplectic form, the introduction 

of (i/>, i/?^) as complex coordinates on the Dirac sector of phase space results in the Poisson brackets 


{F,G}o = -f^ 


dx 


6F 6G 


6G 6F 


(63) 


dipix) (a;) dipix) (a;) J ' 

On a fixed gravitational background, the result of the Legendre transformation would be to drop 
the first term in @ and F[° is the negative of what remains. For the E-D system, we find that 
the final term in (51) contributes to the conjugate momentum of the vierbein, leading to 

pF := -2K\E\ + (64) 

and once again the “geometrodynamical momentum” is not equal to the “gravitational momen¬ 
tum” . However unlike the Maxwell case, the definition of the extrinsic curvature on phase space is 
undisturbed {Kah = and (|7^ ) is unaltered); we find only a contribution to the frame rotation 
generators 

J[ab] ■= -PlaiE\] + \^p^aab'^p■ (65) 

(Note that it is possible, and indeed somewhat simpler, to represent these rotation operators using 
5“ = making use of the spin matrices := = ct“ 1 where (t“ are the Pauli spin 

matrices. We do not do so here simply because we wish to present results that are as similar as 
possible to those of the GL(3,M) case.) In addition to generating infinitesimal frame rotations on 
the GR sector of phase space, we hnd that 


( 66 ) 


(67) 


( 68 ) 


(69) 


J[(h] 

generates frame rotations on the Dirac sector 

and the 0[(3,M) algebra (^) is reduced to that of so(3) 

{J[N],J[M]}= [ 

JY, 

The resulting E-D Hamiltonian is once again in standard form 

H = J d^x {N'rC + > 

where iG = F JF’ and 2ia = -|- 2ia (the GR contributions given by the primed con¬ 

straints (45) of @ with replaced by where 

= + = 77° = VaV-= , (70) 

arise from the Legendre transformation of (^Tj). Note that the equivalence of these with the E- 
D stress-energy tensor components (55) is achieved by making use of the Dirac equation (^) to 
remove time derivatives of the spinor. Also note that the presence of the derivative of the Dirac 
frame rotation generator in the momentum constraint is consistent with the relationship between 
the unprimed and primed constraints in (45b) of [p. From (|70|) we find 

= £^1^], ( 71 ) 

where £' acts as described in p: 

i:L[v>] = E^/^rvam + = E^/^rvam + iv,[/]>. (72) 

Note that we have defined infinitesimal spatial diffeomorphisms to act on the frame as: 2ic[/^]} = 

—Va[/]“£'*b; if one mixes in a spatial frame rotation so that only the symmetric part V 

pears, then the resulting generator acting on ^|J reproduces the action that is used, for example. 
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From (^) one finds Hamilton’s equations for ijj: 

- iVa“)V,[i/;] + i(VjA^]“ - a“Va[iV])^A - (73) 

which are equivalent to (|^). Again making use of the definitions (0) we find (we have used 
= -2e“'’"75) 


wMf] = -mi - cry'll;] + (74a) 

Y\[f] = (74b) 


and using (55) to find 


Tab = V b)V’ + ^6^acKb)d'il^^cr'"^ilJ + ^Sabm'ip'^P'ip, 


(75) 


we find that contributions to (^) and (|7^ are generated to properly reproduce the E-D field 
equations. The constraint algebra combines (|6^) with (^), replacing —> 6ab everywhere in the 

latter. (Note that the derivation of the constraint algebra in |^] could be generalized to the Einstein- 
Dirac system along the lines of [^]; note however the different parameterization of the spinor part 
of phase space and the resulting difference in the form of (|70|) .) 


V. Discussion 

The intent of the preceding and present paper is to describe the Hamiltonian dynamics 
of general relativistic (vacuum and matter coupled) systems with respect to moving frames. This 
extended system of dynamical spatial metric and frame fields encompasses the following two natural 
limits: If we choose the diffeomorphism constraints to act on the components of tensors with 
the usual Lie action and fix the spatial frame to be a coordinate frame, we recover the standard 
coordinate frame approach to canonical general relativity. If we choose the spatial metric to be equal 
to the unit matrix and mix the diffeomorphism generators with the frame transformation generators 
to be compatible with this choice, we recover the orthonormal frame approach to canonical general 
relativity. Clearly there are also a variety of intermediate choices available, however this enlarged 
arena provides a bridge over which one may pass from the standard coordinate frame approach 
to the Lorentz (or orthonormal) frame approach that is more natural for the description of the 
Einstein-Dirac system. 

In this generalized structure, we find that matter fields naively appear to be derivative-coupled 
since the momenta conjugate to the spatial frame fields are not (in general) independent of the 
matter fields. This sort of derivative-coupling turns out to be benign since the field equations 
are equivalent to the coordinate frame approach to the model in question. Indeed this feature is 
a necessary part of the frame approach, and will appear whenever the matter fields behave non- 
trivially under frame transformations. In this way we understand that the derivative-coupling in 
the Einstein-Dirac system is a necessary feature, merely reflecting the need to work with local 
orthonormal frames in order to define the system. We also end up with a simpler description of the 
Einstein-Dirac system than that given in Reference (note however their use of anti-commuting 
fermions), and a much more natural system than that found in Reference m where primary 
importance is still given to the coordinate components of the metric tensor. A more detailed 
examination of the Einstein-Dirac system in adjoint form (in the spirit of Q) is underway, and by 
adopting harmonic coordinate conditions and choosing the frame to fix the Local Lorentz gauge, 
we are able to prove local existence and uniqueness results. 
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